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Abstract. The 222 000 J-band light curves of variable stars detected by the OGLE-II survey in the direction of the Galactic 
Bulge have been fitted and have also been correlated with the DENIS and 2MASS all-sky release databases and with lists of 
known objects. Lightcurves, and the results of the lightcurve fitting (periods and amplitudes) and DENIS and 2MASS data are 
presented for 2691 objects with J-band semi-amplitude larger than 0.45 magnitude, corresponding to classical Mira variables. 
The Mira period distribution of 6 fields at similar longitude but spanning latitudes from —1.2 to —5.8 are statistically indistin- 
guisable indicating similar populations with initial masses of 1.5-2 Mq (corresponding to ages of 1-3 Gyr). A field at similar 
longitude at b — —0.05 from Glass et al. (2001) does show a significantly different period distribution, indicating the presence 
of a younger population of 2.5-3 Mq and ages below 1 Gyr. The A'-band period-luminosity relation is presented for the whole 
sample, and for sub-fields. The zero point depends on Galactic longitude. Simulations are carried out to show that the observed 
dependence of the zero point with I, and the number of stars per field are naturally explained using the model of disk and bulge 
stars of Binney et al. (1997), for a viewing angle (major-axis Bar - axis perpendicular to the line-of-sight to the Galactic Centre) 
of 43 ± 17 degrees. The simulations also show that biases in the observed zero point are small, < 0.02 mag. A comparison 
is made with similar objects in the Magellanic Clouds, studied in a previous paper. The slope of the PL-relation in the Bulge 
and the MCs agree within the errorbars. Assuming the zero point does not depend on metallicity, a distance modulus differ- 
ence of 3.72 between Bulge and LMC is derived. This implies a LMC DM of 18.21 for an assumed distance to the Galactic 
Centre (GC) of 7.9 kpc, or, assuming a LMC DM of 18.50, a distance to the GC of 9.0 kpc. From the results in Groenewegen 
(2004) it is found for carbon-rich Miras that the PL-relation implies a relative SMC-LMC DM of 0.38, assuming no metallicity 
dependence. This is somewhat smaller than the often quoted value near 0.50. Following theoretical work by Wood (1990) a 
metallicity term of the form A/k ~ /3 log Z is introduced. If a relative SMC-LMC DM of 0.50 is imposed, f3 = 0.4 is required, 
and for that value the distance to the GC becomes 8.6 ± 0.7 kpc (for a LMC DM of 18.50), within the errorbar of the geomet- 
ric determination of 7.9 ± 0.4 kpc (Eisenhauer et al. 2003). An independent estimate using the absolute calibration of Feast 
(2004b) leads to a distance estimate to the GC of 8.8 ± 0.4 kpc. 
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1. Introduction 

In the course of the micro lensing surveys in the 1990's the 
monitoring of the Small and Large Magellanic Clouds has re- 
vealed an amazing number and variety of variable stars. A 
big impact was felt and is being felt in any area of variable 
star research, like Cepheids and RR Lyrae stars. Also in the 
area of variability in red variables (RVs) and AGB stars there 
has been remarkable progress. Wood et al. (1999) and Wood 
(2000) were the first to identify and label different sequences 
"ABC" thought to represent the classical Mira sequence ("C") 
and overtone pulsators ("A,B"), and sequence "D" which is not 
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* Table 1 is available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ Figure 1 and the 
Appendices are available in the on-line edition of A&A. 



yet satisfactorily explained (Olivier & Wood 2003, Wood et al. 
2004. Stars on these sequence are referred to as having Long 
Secondary Periods-LSPs). This view has subsequently been 
confirmed and expanded upon by Noda et al. (2002), Lebzelter 
et al. (2002), Cioni et al. (2003), Ita et al. (2004) and Kiss 
& Bedding (2003, 2004), Fraser et al. (2005). These works 
differ in the source of the variability data (MACHO, OGLE, 
EROS, MOA), area (SMC or LMC), associated infrared data 
(Siding Spring 2.3m, DENIS, 2MASS, SIRIUS), and selection 
on pulsation amplitude or infrared colours. In a recent paper, 
Groenewegen (2004; hereafter G04) analysed the OGLE data 
in the SMC and LMC, and correlated these sources with the 
DENIS and 2MASS surveys. The paper discussed the variability 
properties of three samples: about 2300 spectroscopically con- 
firmed AGB stars, around 400 previously known LPV variables 
and about 570 candidate dust-obscured AGB stars. 

The present paper is an extension of the analysis in G04 to 
the OGLE data in the direction of the Galactic Bulge (GB). Also 
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for this area of the sky, several papers exist that use the results 
of the micro lensing surveys and have extended previous clas- 
sical works on Bulge variable stars, like those of Lloyd Evans 
(1976), Glass & Feast (1982), Whitelock et al. (1991), Glass et 
al. (1995; hereafter GWCF), Alard et al. (1996) and Glass et al. 
(2001). 

Alard et al. (2001; herafter ABC01) correlated ISOGAL 
sources within the NGC 6522 and Sgr I Baade windows with 
the MACHO database and present a list of 332 stars with com- 
plete 4-band V, R and [7], [15] magnitudes. Schultheis & Glass 
(2001) extended Alard et al. by also considering the DENIS and 
2MASS data in those fields in general, and for the variables in 
particular. Glass & Schultheis (2002; hereafter GS02) investi- 
gated a sample of 174 M-giants in the NGC 6522 Baade win- 
dow and correlated them with DENIS ISOGAL and MACHO. 
Many stars of spectral type M5 and all M6 and later show vari- 
ation, whereas subtypes M1-M4 do not (see also Glass et al. 
1999). 

Glass & Schultheis (2003; hereafter GS03) investigated 
the variable stars in the NGC 6522 Baade's window using 
MACHO data, and also used DENIS IR data. Of the 1661 se- 
lected stars 1085 were found to be variable. They present K- 
band Pi-relations for sequences "ABCD". 

Wray et al. (2004) investigated small amplitude red giants 
variables in a sub-set of 33 OGLE fields. They identified two 
groups that seem to correspond to groups "A-" and "B-" in Ita 
et al. (2004; also see G04). 

In our paper we describe results on Mira variables selected 
from OGLE Bulge fields. The paper is structured as follows. 
In Section 2 the OGLE, 2MASS and DENIS surveys are de- 
scribed. In Section 3 the model for the lightcurve analysis is 
briefly presented. In the remaining of the paper different results 
are described. The Period-Luminosity diagram is discussed in 
Section 5. A description of the Mira population in respect to the 
overall bulge population is given in Section 8. In Section 9 we 
show that the Miras are distributed in a bar-like structure and 
give the orientation. In the final section we give the distance to 
the Galactic Centre, based on the period-luminosity relation. 

2. The data sets 

The OGLE-II micro lensing experiment observed fourty-nine 
fields in the direction of the GB. Each field has a size 14.2' x 57' 
and was observed in BVI, with an absolute photometric accu- 
racy of 0.01-0.02 mag (Udalski et al. 2002). TableHlists the 
galactic coordinates of the field centers and the total number of 
sources detected in these fields. 

Wozniak et al. (2002) present a catalog of about 222 000 
variable objects based on the OGLE observations covering 
1997-1999, applying the difference image analysis (DIA) tech- 
nique on the /-band data. The data files containing the I- 
band data of the candidate variable stars was downloaded 
from the OGLE homepage (http://sirius.astrouw.edu.pl/~ogle/ 1. 
According to Wozniak et al., the level of contamination by spu- 
rious detections is about 10%, but we presume this level is 
much less at the brighter magnitudes of the LPVs considered 
here. Table|4]lists the number of detected variable stars per field 
(Wozniak et al. 2002). 



Table 2. First entries in the electronically available table, 
which list: OGLE-field and number, other names, spectral 
type and references. ISOGAL sources from the official catalog 
(OGA03) have the prefix "ISOGAL". ISOGAL sources studied 
in OOS03 have the prefix "OOS03". References in column 4 
are given in the bibliography of the main text. 



OGLE-name Other Name 



Chemical type Remarks 



bul.scOI.0098 
bul.sc01.0100 
bul.sc01.0103 
bul.sc01.0128 
bul.sc01.0177 
bul.sc01.0193 
bul.sc01.0235 
bul.sc01.0237 
bul.sc01.0388 
bul.sc01.0558 
buLscO 1.0602 
bul.sc01.0611 
bul.sc01.0878 
bul.sc01.0916 
bul.sc01.1241 
bul.sc01.1253 
bul.sc01.1254 
bul.sc01.1467 
bul.sc01.1533 
bul.sc01.1561 
bul.sc01.1616 
bul.sc01.1738 
buLscO 1.1818 



V4713 Sgr 



BW6_V1_MISC 



BW5.V1.MISC 



BW5.V4.MISC 



The DENIS survey is a survey of the southern hemisphere in 
IJK S (Epchtein et al. 1999). The second data release available 
through ViZier was used (The DENIS consortium, 2003). The 
221801 OGLE objects were correlated on position using a 3" 
search radius and 59894 matches were found. 

The 2MASS survey is an all-sky survey in the JHK S near- 
infrared bands. On March 25, 2003 the 2MASS team released 
the all-sky point source catalog (Cutri et al. 2003). The easiest 
way to check if a star is included in the 2MASS database is by 
uplinking a source table with coordinates to the 2MASS home- 
page. Such a table was prepared for the 221801 OGLE objects 
and correlated on position using a 3" radius. Data on 182361 
objects were returned. 

3. Lightcurve analysis 

The model to analyse the lightcurves is described in detail in 
Appendices A-C in G04. 

Briefly, a first code (see for details Appendix A in G04) se- 
quentially reads in the /-band data for the 222 000 objects, de- 
termines periods through Fourier analysis, fits sine and cosine 
functions to the light curve through linear least-squares fitting 
and makes the final correlation with the pre-prepared DENIS 
and 2MASS source lists. All the relevant output quantities are 
written to file. 

This file is read in by the second code (see for details 
Appendix B in G04). A further selection may be applied (typ- 
ically on period, amplitude and mean /-magnitude), multiple 
entries are filtered out (i.e. objects that appear in different OGLE 
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Table 1. First entries in the electronically available table, which lists: OGLE-field and number, the three fitted periods with errors 

and amplitude (0.00 means no fit), mean I gio, and associated DENIS UK photometry with errors, and associated 2MASS JHK 

photometry with errors (99.9 and 9.99 means no association, or no value). 

OGLE-name Period Amp Period Amp Period Amp J og i e ^denis ^denis -^denis ^2MASS #2MASS ^2MASS 

buLsc01_0098 124.1 ± 0.035 0.322 247.1 ± 0.038 1.520 0.0 ± 0.000 0.000 12.55 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 9.06 ± 0.03 8.06 ± 0.04 7.56 ± 0.02 

buLsc01_0100 332.5 ± 1.306 0.856 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 19.31 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 

bul_sc01.0103 224.2± 0.300 0.082 337.9 ± 0.055 0.904 0.0 ± 0.000 0.000 11.73 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 8.86 ± 0.03 7.80± 0.03 7.11 ±0.02 

bul_sc01.0128 324.0 ± 2.263 0.708 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 19.56 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 

buLsc01_0177 332.3 ± 0.078 2.385 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 14.17 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 8.82 ± 0.03 7.59 ± 0.03 6.90 ± 0.03 

buLsc01_0193 341.1 ± 0.110 0.769 267.1 ± 0.605 0.070 0.0 ± 0.000 0.000 12.60 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 9.12 ± 0.02 7.93 ± 0.02 7.21 ±0.02 

buLscOl-0235 98.4 ± 0.050 0.070 204.7 ± 0.615 0.114 261.8 ± 0.123 0.884 11.79 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 9.17 ± 0.03 8.18 ± 0.05 7.58 ±0.03 

buLscOl-0237 222.5 ± 2.146 0.345 275.0 ± 1.600 0.721 0.0 ± 0.000 0.000 18.66 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 

bul_sc01.0388 107.2 ± 0.192 0.629 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 19.49 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 

bul_sc01.0558 115.4± 0.030 0.167 232.2 ± 0.022 1.186 0.0 ± 0.000 0.000 13.66 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 10.12±0.02 8.92 ± 0.05 8.08 ± 0.02 

bul_sc01.0602 105.6± 0.062 0.114 128.9 ± 0.102 0.106 315.7 ± 0.048 1.498 12.32 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 8.48 ± 0.02 7.43 ± 0.05 6.86±0.02 

bul_sc01.0611 305.8 ± 1.143 0.915 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 19.37 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00±9.99 99.00 ± 9.99 99.00 ± 9.99 

buLscOl-0878 129.1 ± 0.023 0.478 155.4 ± 0.204 0.072 471.0 ± 0.974 0.094 12.34 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 10.06 ±0.03 9.16 ± 0.03 8.73 ± 0.02 

buljsc01_0916 375.4 ± 0.057 0.686 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 11.41 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 7.20 ± 0.02 5.97 ± 0.03 5.28 ±0.02 

bul_sc01.1241 280.1 ± 0.034 0.799 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 11.79 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 8.94 ± 0.04 7.91 ± 0.05 7.31 ±0.02 

bul_sc01.1253 291.9 ± 1.259 0.746 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 19.03 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 

bul_sc01.1254 295.2 ± 1.925 0.995 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 19.81 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 

bul_sc01.1467 129.9 ± 0.042 0.205 167.1 ± 0.063 0.234 519.3 ± 0.106 1.125 12.34 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 8.40 ± 0.02 7.29± 0.06 6.52±0.02 

bul_sc01.1533 116.2± 0.008 1.156 336.2± 0.360 0.181 505.8 ± 0.855 0.148 12.85 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 10.39 ±0.04 9.46 ± 0.05 8.97 ±0.03 

bul_sc01_1561 46.8 ± 0.056 0.558 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 18.22 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 

buLsc01_1616 457.7 ± 0.104 0.903 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 12.33 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 8.81 ± 0.02 7.25 ± 0.04 6.38 ±0.02 

bul_sc01.1738 155.0± 0.055 0.153 458.0 ± 0.267 1.143 0.0 ± 0.000 0.000 13.27 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 9.85 ± 0.02 8.31 ± 0.07 7.32±0.02 

bul_sc01.1818 279.3 ±0.052 0.785 0.0 ± 0.000 0.000 0.0 ± 0.000 0.000 11.61 99.00 ± 9.99 99.00 ± 9.99 99.00 ± 9.99 9.22 ± 0.02 8.21 ± 0.05 7.56±0.03 



fields), and a correlation is made with pre-prepared lists of 
known non-LPVs and known LPVs or AGB stars. The output 
of the second code is a list with LPV candidates. 

The third step (for details see Appendix C in G04) consists 
of a visual inspection of the fits to the light curves of the candi- 
date LPVs and a literature study through a correlation with the 
SIMBAD database. Non-LPVs are removed, and sometimes the 
fitting is redone. The final list of LPV candidates is compiled. 

Details on the small changes in the codes w.r.t. the imple- 
mentation in G04 are given in Appendix A of the present paper. 



4. Comparison of the datasets 

4. 1. Astrometry 

As in G04, the correlation between the OGLE objects and 
known LPVs and AGB stars, and known non-LPVs, is actu- 
ally done in 2 steps. In the first step the correlation is made (for 
a 3" search radius), and the differences and spread in ARA 
cos((5) and A5 are determined. These mean offsets are then ap- 
plied in most cases to make the final cross-correlation, and this 
usually increases the number of matches. The results are listed 
in TablefJ 

4.2. Photometry 

As in G04, a comparison was made between the (mean) OGLE 
/ and the (single-epoch) DENIS /, and between the (single- 
epoch) DENIS JK and the (single-epoch) 2MASS JK mag- 
nitudes. This was done by selecting those objects with an am- 
plitude in the /-band of < 0.05 mag. 

Figure |2] shows the final results when offsets /(denis- 
ogle) = —0.03, J(denis-2mass)= —0.02, and A"(denis-2mass)= 
—0.03 are applied. The latter values are consistent with those 
derived in OOS03 based on the 2MASS second incremental 
data release who found J(denis-2mass)= —0.02 ± 0.09, and 
J fY(denis-2mass)= -0.00 ± 0.07. 




14 
I (ogle) 




12 13 
J (2MASS) 




■ 



K (2MASS) 

Fig. 2. Difference in photometry, after the following offsets 
have been applied: /(denis-ogle) = —0.03, J(denis-2mass)= 
-0.02, if(denis-2mass)= -0.03. 



5. Period-Luminosity relations 

The full machinery outlined in Section 3 was performed. As in 
G04, all derived periods are given in Tableland are shown in 
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bul_sc01_0103 




Fig. 1. First entries of electronically available figure with all lightcurves. The fit is indicated by the (red) solid line. Crosses 
indicate data points not included in the fit. 



Figure^ The present paper discusses only objects which have 
at least one period with an /-band amplitude larger than 0.45 
magnitudes 1 , i.e. classical Mira variables (e.g. Hughes 1989). 
After visual inspection of the lightcurves a sample of 2691 such 
objects remain. The number of objects per field is listed in the 
last column of Tab. [4] 

Table^lists the stars in the sample, the fitted periods with 
errors and amplitudes, and the DENIS and 2MASS photometry 
of the associated sources. Table |2] lists alternative names and 
references from the literature. Figure^presents lightcurves and 
their fits. 

In the discussion that follows, magnitudes are de-reddened 
using the Ay values that correspond to the respective OGLE 



1 The amplitude, A, is used in the mathematical sense in the present 
paper, y = A sin x. The peak-to-peak amplitude is 0.90 mag. 



field taken from Sumi (2004; and Ay = 6.0 for the field 44 that 
they do not discuss), and selective reddenings of Aj/Ay = 
0.49, Aj/A v = 0.27, A n /Ay = 0.20, A K /Ay = 0.12 
(Draine 2003) and implicitly assuming that all objects suffer 
this reddening value (i.e. ignoring differential reddening within 
a field, and ignoring that foreground and background objects 
would suffer a different reddening). Sumi's method is based 
on the absolute magnitude of the Red Clump giants and is ab- 
solute calibrated using the (V — K) colours of 20 RR Lyrae 
stars in Baade's window. Popowski et al. (2003) present an ex- 
tinction map (over 9000 resolution elements of 4x4 arc minute 
size) towards the GB based on MACHO V, R photometry, un- 
der the assumption that colour-magnitude diagrams would look 
similar in the absence of extinction. For the centre of the OGLE 
fields it was checked if there was a tile in the Popowski et al. set 
within 0.05 degrees distance. For those, the value of the visual 
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Table 3. Comparison of coordinates, and number of positional matches, before and after a correction was applied. 



OGLE 

-(other) 


ARA cos(<5) 


A<5 


N 


ARA cos(<5) 


A<5 


N 


remark 


IRAS 


0.14± 0.94 


0.50± 1.44 


80 








IRAS sources, not corrected 


OGA03 


-0.80± 0.85 


-0.16± 0.89 


1309 


-0.11± 0.82 


-0.02± 0.90 


1319 


ISOGAL sources from Omont et al. (2003; OGA03) 


OOS03 


-0.80± 0.88 


-0.30± 0.95 


800 


-0.10± 0.86 


-0.03± 0.99 


817 


ISOGAL sources from Ojha et al. (2003; OOS03) 


GS02 


-0.02± 0.58 


0.30± 0.61 


79 


-0.04± 0.67 


0.00± 0.61 


80 


MACHO sources from Glass & Schultheis (2002; GS02) 


ABC01 


-0.82± 0.98 


-0.70± 0.96 


293 


-0.17± 1.08 


-0.12± 1.01 


320 


ISOGAL sources from Alard et al. (2001; ABC01) 


OGLE-I 


-0.08± 0.70 


-0.18± 0.44 


728 








OGLE-I sources, not corrected 


BMB 


-0.04± 0.70 


-0.52± 0.90 


282 


+0.10± 0.79 


-0.08± 0.91 


286 


sources from Blanco et al. (1984; BMB) 


TLE 


-0.08± 0.84 


-0.25± 1.43 


21 








sources from Lloyd-Evans (1976), not corrected 


B84 


-0.46± 0.75 


0.50± 0.97 


61 


-0.02± 0.76 


0.08± 1.05 


63 


sources from Blanco (1984; B84) 



extinction has been listed, next to the value in Sumi in Tabled 
The rms difference Ay for the 21 field with values from both 
references is 0.18. Finally, Schultheis et al. (1999b) presented 
a reddening map for the inner GB comparing DENIS J, K pho- 
tometry to isochrones. Table|4]lists for two OGLE fields the val- 
ues they find: SC44 which was not considered by Sumi (2004), 
and SC5 for which Sumi derive a larger Ay than Schultheis et 
al.: 5.73 versus 4.13. 

In the further discussion we only use periods that fulfil the 
following conditions are used in the calculations (with AP the 
error in the period): AP/P < 0.01 for P < 500 d ; AP < 5 d 
for 500 d < P < 800 d and AP < 1.5 d for P > 800 d . The 
latter constraint is necessary because the long periods become 
comparable to the length of the dataset. 

Figure [5] shows the LT-band PL-relation for all periods 
which have an /-band amplitude larger than 0.45 magnitude 
and (J - K) < 2.0 among the 2691 stars. The cut in (J - K) 
colour is needed to prevent that the LT-magnitude is affected 
by circumstellar extinction, as shown in G04. Like G04, the 
LT-magnitude is on the 2MASS system, and is the average of 
the DENIS and 2MASS photometry. In particular, if both DENIS 
and 2MASS LT-band data is available, the DENIS data point is 
corrected as explained above (i.e. 0.03 mag added), and aver- 
aged with the 2MASS data point. This should take out some of 
the scatter in the Pi-diagram, as the effect of the variability in 
the L'-band is reduced. If only DENIS is available, the corrected 
value is used. In the left-hand panel the boundaries of the boxes 
"A-, A+, B-, B+, C, D" have been taken from G04, but shifted 
by —4.0 to account for the approximate difference in distance 
modulus (DM), as, e.g., follows from the recent determination 
of 7.9 ± 0.4 kpc (Eisenhauer et al. 2003) for the distance to the 
GC, and 18.50 for the DM to the LMC (e.g. recent reviews by 
Walker 2003, Feast 2004a). 

There is a reasonably well defined sequence in Box "C", 
but when compared to the similar figure for the SMC and LMC 
in G04 (his figure 3) some differences can be remarked as well. 
In particular, for the present Bulge sample there are a few ob- 
jects located in Box "B+", and in particular many in Box "D". 
In the SMC and LMC, for this cut in amplitude, there are none 
in Box "B+" and few in "D". Several issues may play a role. 
Applying a certain cut in amplitude may sample slightly differ- 
ent variables in SMC, LMC and Bulge. Figure 3 in G04 clearly 
shows how lowering the cut in amplitude results in a populating 
Box "B+" and then "A+", and increases the number of objects 



in "D". Another effect is the possible contribution of objects 
in the foreground and background of the Bulge, the depth of 
the Bulge, and fourthly, the orientation of the Bar, as the OGLE 
fields span 20 degrees in longitude (this last effect will be dis- 
cussed later). Finally, the difference in DM may be different 
from the adopted value of 4.0. 

To verify if the objects in Box "D" actually show LSP, they 
were all visually inspected, and in fact few have, in agreement 
with the finding for LMC and SMC (for amplitudes >0.45 
mag). This would call for a enlargement of Box "C" to properly 
sample the PL-relation of the large amplitude (Mira) variables. 
To define this enlarged box the PL-relation was inspected for 
each field independently. The right-hand panel in Figure [5] 
shows the finally adopted boundaries of Box "C", which im- 
plies that Box "D" has contracted. Stars inside this redefined 
Box will be used to define the PL-relation. The LT-band PL- 
relation is determined to be: 

m K = (-3.37 ± 0.09) log P + (15.44 ± 0.21) (1) 

with an rms of 0.42 and based on 1292 stars, and is shown in 
Fig. [5] The value of the slope is consistent with the median 
value when the PL-relation is determined for all fields individ- 
ually. 

For reference, fitting all stars in Figure [3] for a fixed slope 
of -3.37 results in a ZP of 15.47 ± 0.55. 

6. Historical versus current periods 

Table |3 compares the period derived in the present paper (the 
one with the largest amplitude) with values derived in the liter- 
ature. There are three cases where a previously determined pe- 
riod may be a harmonic of the present period but overall there 
is good agreement between periods. In the 12 cases where there 
is a period available from LE76 (with the photographic plates 
taken between 1969 and 1971, hence 28 years of time differ- 
ence with OGLE) there is no clear case for a star that changed 
period. By comparison, G04 found that about 8% of LMC vari- 
ables changed their period by more than 10% over a about 17 
year timespan. To find out of 12 in the present sample is con- 
sistent with this. 

7. Colour-colour diagrams 

The 2MASS and DENIS Colour-colour and colour-magnitude 
diagrams are shown in Fig. |4] together with that of spectro- 
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Fig. 3. If -band Pi-relation for periods with an /-band amplitude larger than 0.45 mag and (J — K)o < 2.0. Left panel: 
Sequences/Boxes "A+, A-, B+, B-, C, D" taken from G04 for the LMC and shifted by -4.0 in distance modulus. Right panel: 
enlarged Box "C" (and contracted Box "D",) to define the region for which the PL-relations will be computed. Only periods 
from TableQ]that fulfil AP/P < 0.01 for P < 500 d ; AP < 5 d for 500 d < P < 800 d and AP < 1.5 d for P > 800 d are plotted 
and used in the analysis, known M-stars are indicated by open circels. The line in the right-hand panel indicates the PL-relation 
of Equation 1 . 



scopically confirmed M-stars in the LMC (see also Figure 12 
in G04). There appear to be more redder stars in the Bulge 
sample, but this is likely due to a under representation in the 
LMC sample as this was restricted to spectroscopically known 
M-stars (i.e. in general optically selected). The sample of can- 
didate infrared-selected AGB stars in the LMC [Figure 12 in 
G04] does cover the (/ — K) and (J — K) colour range ob- 
served in the Bulge). 

The other main difference is that the Bulge stars are redder 
by ~ 2 mag in (J— J)q compared to both LMC and SMC stars, 
as was also shown by Lebzelter et al. (2002) in a comparison 
of LMC and Bulge variable stars. As the diagrams involving 
J, H, K colours appear similar, it seems that this difference in 
(I — J) must be due to a difference in I. The /-band mea- 
surements of M stars is strongly affected by the TiO and VO 
molecular absorption features (Lancon & Wood 2000). It is ex- 
pected that for larger metallicities these lines will be stronger 
(Schultheis et al. 1999a) which will lead to redder (I - J) 
colours. 

The bullets connected by a line in the Bulge DENIS (/ — 
J) — (J — K) colour-colour diagram are the average colours of 
Ml, M2, .., M6, M6.5, M7, M8 giants in the NGC 6522 Baade's 
window (Blanco 1986, GS02). There is a spread of typically 
0.3-0.5 mag in (/ — J) and 0.2-0.3 mag in (J — K) around 
these means, and there is only 1 M8 giant in their sample. The 
colours of the Miras follow those of normal giants well until 
M6.5, when the Miras become redder in (I — J). 

There are also stars redder in (J — K) than the single M8 
star in the sample of GS02, indicating either the presence of 
later spectral types or the on-set of circumstellar reddening. 

The conclusion of GS02 that "Many M5 and all stars M6 
and later show variation, whereas subtypes (M1-M4) do not" is 
confirmed here, as there are essentially no objects located in the 
region of the DENIS (I — J) — (J — K) colour-colour diagram 
occupied by spectral types of M4 and earlier. 



8. Mira bulge population as function of latitude 

Figure [5] shows the period distribution of Miras in Box "C". 
A distinction is made between all Miras and those with (J — 
K)o < 2.0 (dashed histograms). The latter selection minimises 
any influence of circumstellar extinction. For comparison, the 
period distribution of LMC and SMC Miras is also shown 2 . 
The Kolmogorov-Smirnov (KS) test is performed to indicate 
that the probability that the period distributions of Bulge-LMC, 
Bulge-SMC, LMC-SMC are the same for all stars (those with 
(J - K) < 2, respectively) is, respectively 0.36 (10 -8 ), 0.05 
(0.31) and 0.05 (0.05). 

Any difference, in particular between Bulge and LMC pe- 
riod distribution, is difficult to quantify further as this de- 
pends in a complicated way on the Star Formation History 
and evolutionary tracks (T e g - Luminosity - Mass - metallic- 
ity). Regarding the period distribution of Bulge Miras as such, 
previous studies are limited to selected small fields (e.g. TLE, 
GWCF, Glass et al. 2001). Whitelock et al. (1991) present the 
period distribution of about 140 IRAS sources but no direct 
comparison is possible because of the difference in the selec- 
tion criteria of the samples. 

Figure |6] shows the period distribution of selected fields 
with very similar longitudes that cover a range in latitudes (the 
stars with (J — K)q < 2.0 are shown as dashed histogram 
again). To add a field even closer to the GC than surveyed by 
OGLE the data in Glass et al. (2001, 2002) is considered on a 
field centered on / = —0.05, b = —0.05. They present the re- 
sults of a LT-band survey of 24 x 24 arcmin 2 for LPVs down 
to K ~ 12.0. From the list of 409 stars, 14 were removed be- 
cause of double entries, quality index of zero, uncertain or no 
listed period or amplitude. The coordinates were uploaded to 
the IPAC webserver and 2MASS data within 2.5" was retrieved 

2 Derived following the implementation of the code and definition 
of the "boxes" as in G04, and applying the same selection criteria as 
in the present paper, i.e. /-amplitude larger than 0.45 mag. 
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Fig. 4. Colour-Magnitude and colour-colour diagrams using 2MASS and DENIS photometry for the Bulge stars (left), and spectro- 
scopically confirmed M-stars in the LMC (right, fromG04). The bullets connected by aline in the Bulge DENIS (7 — J) — (J — K) 
diagram are the average colours of Ml, M2, .., M6, M6.5, M7, M8 giants in the NGC 6522 Baade's window (Blanco 1986, Glass 
& Schultheis 2002). 



for these 395 stars to get information on the (J — K) colour. 
As an additional check, and to eliminate multiple stars within 
the search circle, it was verified that the single-epoch 2MASS 
K magnitude is consistent with the mean if -magnitude and 
amplitude listed in Glass et al. For 345 stars 2MASS data is 
available. The magnitudes are corrected for interstellar redden- 
ing using the extinction map of the inner GB at 2' resolution by 
Schultheis et al. (1999b). The extinction value of the nearest 
available grid point in this map is taken. The extinction val- 
ues range between 18.5 and 30.4 with a mean of 24.7. The top 
panel in Figure [6] list 333 stars with K > 7 (to eliminate 3 
very likely foreground objects), and if-band amplitude larger 
than 0.35 (to correspond roughly to the cut in /-band amplitude 
of 0.45 mag), and 88 (the histogram with slanted hatching), or 
236 (dashed histogram) stars which also have (J — K )o < 2.0. 
The latter sample is the one that results when the reddening 
values from Schultheis et al. are multiplied by 1.35. They men- 
tion themselves that the reddening may be underestimated in 
the direction of the GC because of J-band non-detections. For 
the one field in common, their value is a factor 1.3-1.4 smaller 
than derived by Sumi (2004). In addition, for their default red- 
dening (the histogram with slanted hatching in Figure |6j there 
would be many stars even at periods shorter than about 250 
days which still would have ( J — K)o > 2.0 which is not ob- 
served in the other fields. This could off course be real, but it is 
generally believed (e.g. Launhardt et al. 2002) that the popula- 
tion of low- and intermediate mass stars in the Nuclear Bulge 
(the inner about 30 pc from the GC) and GB are similar, but 
that in the former there is an overabundance of 10 7 — 10 8 year 
old stars. In this picture one would expect the period distri- 
butions to be similar at shorter periods, essentially indepen- 
dent of latitude. Therefore the period distribution of stars with 
(J — K) < 2.0 for the increased reddening is adopted. 



The KS test is performed on consecutive fields in latitude 
for the distributions based on the stars with (J — K)q < 2.0. 
It is found that the probability that the distributions are the 
same is 10~ 10 for the b = -0.05/ -1.21 fields, 0.50 for the 
b = -1.21/ -1.39 fields, 0.80 for the b = -1.39/ -1.81, 
and > 0.99 for the fields at more negative latitudes. The con- 
clusion is that the period distributions of the fields at and be- 
low — 1.2 degree are statistically indistinguisable, but that the 
field at —0.05 latitude has a significantly different period dis- 
tribution (the probability that this distribution is the same as 
the distribution of the combined 6 OGLE fields is 10~ 22 ). This 
conclusion is independent of the assumed reddening of the in- 
ner Bulge field, which influences how many stars will have 
(J — K)q < 2.0. For the default reddening of Schultheis et 
al., the probability that the distributions are the same for the 
b = -0.05/ -1.21 fields is still only 0.0033. The difference in 
the period distributions is especially clear at longer periods. Of 
the 236 stars in the inner field with (J — K)o < 2.0 61 have 
P > 500 days, while in the other fields this is 3 out of 367. 

The difference in period distribution might be due to an 
under representation of short period stars in the inner field. 
However, Figure 4 in Glass et al. illustrates that the expected 
if -magnitudes at short periods are not fainter than the com- 
pleteness limit of their survey. In fact, Glass et al. mention that 
they expect that the number of short-period Miras (P < 250 
days) is at least 75% complete. As a test, one-third of stars 
with P < 250 days were randomly duplicated and added to 
the sample, and the KS test repeated to find again a large dif- 
ference between the period distribution of the field at —0.05 
degrees and the other fields. 

The difference is emphasised in Figure0where the scaled 
period distribution of stars which have (J — K)o < 2.0 in the 5 
fields between b = —1.39° and b = —5.8° has been subtracted 
from the inner field. The scaling was done in such a way that 
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at shorter periods the two distributions would cancel at a level 
of la (based on Poisson errors). Even if the scaling is done in 
a slightly different way the result is always very similar, in the 
sense that there is a significant (> 4cr) overabundance of LPVs 
in the inner field bewteen about 350 and 600 days. 

The conclusion is that there is a significant population of 
LPVs with period >,500 days present in the inner field, which 
remains barely present at latitude —1.2°, and is absent for 
b< — 1.4°. This was indirectly noted by Glass et al. who noticed 
that the average period of the stars in this field at b = —0.05° 
is 427d (and that of the known OH/IR stars 524tf), while the 
average period in the Sgr I window (b = —2.6°) is 333d, with 
no known OH/IR stars (GWCF). 

To quantify the nature of the Mira Bulge population, syn- 
thetic AGB evolutionary models have been calculated, which 
are described in detail in Appendixlcl 

In brief, the synthetic AGB code of Wagenhuber & 
Groenewegen (1998) was finetuned to reproduce the models 
of Vassiliadis & Wood (1993) and then extended to more ini- 
tial masses and including mass loss on the RGB. For several 
initial masses the fundamental mode period distribution was 
calculated for stars inside the observed instability strip and 
when the mass loss was below a critical value to simulate the 
fact that they should be optically visible. Vassiliadis & Wood 
(1993) provide calculations for 4 different metal abundances: 
Z = 0.016,0.008,0.004 and 0.001. We used the models for 
Z = 0.016, representing a solar mix, which are most appro- 
priate for our Bulge sample (e.g. Rich 1998). We also show 
that our results are essentially unchanged if Z = 0.01 or 0.02 
are adopted. From the comparison of the observed period dis- 
tribution for fields more than 1.2° away from the galactic cen- 
tre with the theoretical ones, we deduce that the periods can 
be explained with a population of stars with Main Sequence 
masses in the range of 1.5 to 2.0 M Q . A possible extension to 
smaller masses is possible, but not necessary to explain the pe- 
riods below 200 days. To explain the excess periods in the range 
of 350-600 days observed closer to the centre we need initial 
masses in the range 2.5 - 3 M Q . The presence of more mas- 
sive stars in the inner field at b = —0.05° cannot be excluded, 
as it turns out that for more massive stars the optically visi- 
ble Mira phase is essentially absent. Sevenster (1999) analysed 
OH/IR stars (which are LPVs with longer periods and higher 
mass loss rates than the Miras) in the inner Galaxy and came to 
the conclusion that OH/IR stars in the bulge have a minimum 
intial mass of about 1.3 M Q , based on an analysis of infrared 
colours, compatible with our results. We briefly mention here 
the result from Olivier et al. (2001) who studied a sample of 
LPVs in the solar neighbourhood with periods in the 300 to 
800 days range. They conclude that majority of these stars had 
initial masses in the range of 1 - 2 M , with an average value of 
1 .3 M Q , lower than what we find for the 300 to 600 days range 
sample. This difference may be explained by the fact that our 
conclusions are only valid for a sample with no or only low 
mass loss rates (<,5 x 10 -6 M Q jyr ), contrary to their sam- 
ple which was selected to contain stars with significant mass 
loss (~ 10 -5 M Q jyr). As can be seen in the Vassiliadis & 
Wood (1993) models, the period increases considerably when 
the stellar mass is reduced by the mass loss process. 



We do not see a variation in the period distributions for the 
higher latitude fields (beyond 1° latitude) and can consider this 
as a homogeneous "bulge" population, which according to the 
Vassiliadis & Wood (1993) model has ages in the range of 1 
to about 3 Gyr. The excess population closer to the Galactic 
Centre is younger than 1 Gyr. According to Launhardt et al. 
(2002), the Nuclear Bulge (approximately the central degree) 
contains besides the bulge population seen at higher latitudes 
also an additional population due to recent star formation closer 
to the galactic centre. Blommaert et al. (1998) find that the ex- 
trapolation of the number density of bulge OH/IR stars towards 
the galactic centre would explain half of the galactic centre 
OH/IR population, but that an additional population, intrinsic 
to the galactic centre, exists which agrees with what we see in 
the distribution here. 

The formation history of the Bulge is still a matter of de- 
bate. In several works like in Kuijken & Rich (2002) and re- 
cently in Zoccali et al. (2004), the bulge is considered to be 
old (> 10 Gyr) and formed on a relatively short timescale 
(< 1 Gyr) (e.g. Ferreras et al. 2003). On basis of the mod- 
elling of colour-magnitude diagrams, Zoccali et al. claim that 
no trace is found of any younger stellar population than 10 Gyr. 
The bulge Miras do not fit in this picture as, according to our 
analysis, they are considerably younger. The field studied by 
Zoccali et al. is centered at (l,b) = (0.277, —6.167) so at a 
slighter higher latitude than our extreme fields (b i=s —5.8°). 
Although we are limited by the small number of Miras detected 
at the highest latitude fields, we do not see a change in period 
distribution for those fields. Zoccali et al. (2004) acknowledge 
the presence of Miras, but consider them as part of the old pop- 
ulation. It is true that Miras are also detected in globular clus- 
ters and thus can be associated with old ages, as is the case for 
a 1 Mq star in the Vassiliadis & Wood (1993) model but these 
stars produce periods shorter than 200 days (Figure C. 1), insuf- 
ficient to explain the period distribution seen in the bulge. The 
periods of Miras in Globular Clusters range from 150 to 300 
days (Frogel & Whitelock 1998) and so only overlap with the 
shorter periods of the bulge Miras. 

Our results agree more with the analysis of the infrared 
ISOGAL survey discussed in van Loon et al. (2003). They con- 
clude that the bulk of the bulge population is old (more than 
7 Gyr) but that a fraction of the stars is of intermediate age (1 
to several Gyr). The Miras in our study can thus be considered 
as the intermediate age population seen in their analysis, van 
Loon et al. also see evidence for an even younger population 
(< 200 Myr), but according to our findings, this would be re- 
stricted to the area close to the Galactic Centre. 

Our discussion on the ages of the Mira stars is based on 
the assumption that they have evolved from single stars. An 
alternative scenario as suggested by Renzini & Greggio (1990), 
would be that the brighter (longer period) Miras could evolve 
from close binaries where the components coalesced to form 
one single star. This could lead to an underestimation of the 
age as the Mira essentially is the product of lower mass and 
thus older stars. This scenario may seem in better agreement 
with the idea that the bulge consists of an old stellar population. 
It however suffers from the same problem as the intermediate 
age population, in the sense that also no clear evidence for Blue 
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Fig. 5. Period distribution of large amplitude variables in box 
"C", for Bulge (top), LMC (middle), and SMC (bottom). For 
the shaded histogram only stars with (J — K)o < 2.0 have 
been included. 



Stragglers (which would be the Main Sequence counterpart of 
the Miras) is found (Kuijken & Rich, 2002). 

If indeed the bulk of the bulge population is old and formed 
quickly and if the Miras are of intermediate age, then our Miras 
must be representatives of a population which was added at a 
later stage and it is unclear how it relates to the overall bulge. 
An interesting scenario suggested in Kormendy & Kennicutt 
(2004) is the one in which a secondary bulge or also called 
pseudo-bulge forms within an old bulge. Such a process would 
be connected to the presence of a "bar" which would add 
"disky" material into the old classical bulge. The Miras are in- 
deed situated in a bar-structure as is discussed in the following 
section. 
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Fig. 7. Period distribution of the stars in the b = —0.05° field 
with (J — K)q < 2.0 minus the suitably scaled period dis- 
tribution of the stars in the b = —1.39 to —5.8° fields with 
(J — K )o < 2.0. The scaling is such that below <,300 days the 
two distributions cancel at the 1 a level. 

9. The orientation of the Bar 

Table [6] lists the zero points (ZPs) of the the if-band PL- 
relation (for a fixed slope of —3.37) for the OGLE fields in- 
dividually. To increase the statistics, some neighbouring fields 
have been added together, as indicated in the first column of 
the table. The galactic coordinates listed are the mean values 
of all individual objects, rather than the mean of the field cen- 
tres. Figure |9] plots these ZPs (with error bars) as a function 
of Galactic longitude. There is a clear correlation; the formal 
weighted fit has a slope of —0.023 ±0.005 (magnitude/degree). 
Restricting the fields to those with longitudes —5 < I < +5 
(reducing the contamination by disk stars, see AppendixlBl the 
fit becomes: 



m K = (-0.0192 ± 0.0087) I + (15.484 ±0.019) 



(2) 



with an rms of 0. 10 and based on 32 fields. 

The interpretation of this correlation is that the Bulge Miras 
are located in the Galactic Bar that has a certain orientation 
towards the observer. A similar correlation was found by Wray 
(2004) who concluded that an appropriately chosen ZP in I 
for the small amplitude OGLE variables in their sample (which 
they identify as to correspond to in Box "A-") correlated with 
Galactic longitude. No estimate for the orientation of the Bar 
was given however. 

In Appendix|B]Monte Carlo simulations are carried out in 
order to quantify two issues: can these observations be used 
to constrain the orientation of the Galactic Bar, and, second, 
given the specific location of the OGLE fields, if there is any 
bias in the derived zero point compared to a fiducial ZP when 
all Miras would be located exactly in the Galactic Centre (GC). 
As described in Appendix[B] for a spatial distribution of Bulge 
and Disk stars following Binney et al. (1997), viewing angles 
tfi of 43 and 79 degrees (see the orientation in Figure |8j result 
in slopes (magnitude versus /, Eq. 2) in agreement with obser- 
vations. However, the model with <f> = 43° gives a much better 
fit to the number of stars per field. The bias in the ZPs is es- 
sentially independent of viewing angles, and for the best fitting 
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Fig. 6. Mira period distribution for 7 fields with similar longitudes but a range in latitudes (as indicated in the top right corner). 
For the field at b ~ -5.8°, OGLE fields 6 and 7 have been combined. For the shaded histograms only stars with (J - K)o < 2.0 
have been included. The field at (—0.05, —0.05) is based on Glass et al. (2001), see main text for details. The histogram with 
slanted hatching is for the reddening by Schultheis et al. (1999b) for stars in this field, the shaded histogram for the adopted 
reddening which is 1 .35 times larger. 
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Fig. 8. Schematic drawing of the orientation of the major-axis 
of the Galactic Bar w.r.t. the Sun and the Galactic Centre. The 
z-axis is directed towards the reader. See AppendixlBl 

model the observed ZP derived for all stars (Eq. 1) is too bright 
by 0.018 mag (± 0.013), while the ZP in Eq. 2 is too bright by 
0.002 (± 0.021) mag. 

The preferred value of <fi — 43° is in agreement with the 
values of about 45° by Whitelock (1992), based on 104 IRAS 
detected Mira variables, and the preferred value of 46° by 
Sevenster et al. (1999), based on an analysis of OH/IR stars 
in the inner Galaxy. 

Other values in the literature are usually much larger, be- 
tween 60 and 80 degrees: Dwek et al. (1995) and Binney et 
al. (1997), based on COBE-DIRBE data, Stanek et al. (1997), 
based on bulge red clump stars, Robin et al. (2003) and Picaud 
& Robin, based on colour-magnitude fitting. Sevenster et al. 
(1999), however, argues that these values are commonly found 
when no velocity data is available, the longitude range is too 
narrow or when low latitudes are excluded. It is also possi- 
ble that these studies are tracing other populations, which may 
be differently distributed than the Miras. Whitelock et al. and 
Sevenster et al. do use populations closely related to the Mira 
stars and find an angle of the bar close to the one we derive. 

10. The distance to the Galactic Centre 

A slope of the LT-band PL-relation of -3.37 ± 0.09 is de- 
rived. There appear not to exist many previous determina- 
tion of this quantity. Recently, GS03 derived a PL-relation in 
NGC 6522 based on 34 MACHO variables with r-amplitude 
>1.5 and DENIS K photometry: m K = -4.6 logP + 18.1. 
No errors or rms were given, as-by their own account-this 




Galactic Longitude 

Fig. 9. Zero point of the LT-band PL-relation as a function of 
Galactic longitude. Galactic latitudes below —4.0 are indicated 
by filled triangles, those larger than —2.6 by filled circles, and 
the remaining by open circles. Error bars are also plotted. The 
lines represent linear least-squares fits to all fields (top panel), 
and those with | Z |< 5° (bottom panel). 

fit was made by eye. Much better agreement is found with 
GWCF. Based on multi-epoch data of 55 stars they found 
m K = (-3.47 ± 0.35) logP + (15.64 ± 0.86) (rms=0.35) 
in the Sgr I field. 

Zero points for the LT-band PL-relation have been derived 
in two ways. First, a direct fit to all stars resulting in (15.44 ± 
0.21), and secondly, determining ZPs per (sub)-field, and fitting 
this as a function of I, resulting in (15.484 ± 0.019). 

Applying the small bias corrections discussed at the end of 
Sect. |9] and averaging over the two estimates, the adopted K- 
band PL-relation for Miras at the GC is: 



m K = -3.37 log P+ (15.47 ±0.03) 



(3) 



The derived PL-relation can be compared to the one de- 
rived for 83 O-rich LPVs in the LMC derived in G04: m K = 
(-3.52 ± 0.16) log P + (19.56 ± 0.38), with an rms of 0.26. 
Since the slopes are not exactly the same, the magnitudes are 
compared at the approximate mean period of logP = 2.45. 
The difference in magnitude is 3.72. Adopting the LMC based 
slope of —3.52 for the GB Miras, and re-fitting the ZP, the bias 
corrected ZP would become 15.85, resulting in a GB-LMC DM 
difference of 3.71, essentially the same value. If the distance to 
the GC is assumed to be 7.94 kpc (Eisenhauer et al. 2003; in 
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a recent preprint this was even lowered to 7.62 ± 0.32 kpc, 
Eisenhauer et al. (2005)), then the LMC would be at a DM = 
1 8 .2 1 , or if the DM to the LMC is assumed to be 1 8 .50 (Walker 
2003, Feast 2004a), then the GC would be at 9.0 kpc. A sim- 
ilar result was found by GWCF who derived a distance to the 
GC of 8.9 ± 0.7 kpc, assuming 18.55 for the LMC DM and 
<p = 45°. The analysis so far has assumed no metallicity de- 
pendence of the Mira PL-relation. Wood (1990) present linear 
non-adiabatic pulsation calculations that suggest a dependence 
of the form logP ~ 0.46 log Z + 1.59 log L, but he notes that 
in the if -band the dependence is expected to be weaker and 
following the example he presents one infers a dependence of 
0.25 log Z in the if -band. In G04 if -band PL-relations were 
derived for carbon-miras in the SMC and LMC. At a character- 
istic period of logP = 2.45 one infers a relative difference in 
DM of 0.38, which is smaller than the commonly quoted value 
of near 0.50 (0.48-0.53 ± 0.1 1, FO cepheids [Bono et al. 2002], 
0.46-0.51 ± 0.15, FU cepheids [Groenewegen 2000], 0.44 ± 
0.05, TRGB [Cioni et al. 2000]). This may hint at a metallic- 
ity dependence of the Mira if -band PL-relation. To test this 
hypothesis, a correction to the if -magnitude of +/3\ogZ will 
be assumed 3 (for both O- and C-rich LPVs), and the Bulge, 
LMC, and SMC will be assumed to have solar, solar/2 and so- 
lar/4 metallicity, respectively. For a value (3 — 0.25 the relative 
SMC-LMC DM based on the C-Miras is increased from 0.38 
to 0.46, while the relative DM LMC-GC is increased from 3.72 
to 3.80. If the relative SMC-LMC DM is fixed at 0.50, then 
P = 0.40 is required, and the relative DM LMC-GC becomes 
3.84 for that value. For a LMC DM of 18.50, the distance to 
the GC then becomes 8.6 kpc. The error in this value is some- 
what difficult to estimate as the PL-relations derived in G04 
and here are from-at best-the average of two if values. Work 
by Feast et al. (1989) indicates that in the case of multi-epoch 
data (and for the small depth effect in the LMC) the intrin- 
sic dispersion in the PL-relation is about 0.13 mag. Therefore 
we assign an error of 0.18 to the difference in DM. This im- 
plies an error of 0.7 kpc. Based on this large sample of Mira 
variables in the direction of the GB the conclusion is that the 
distance to the GC is between 8.6 and 9.0 (±0.7) kpc, depend- 
ing on the metallicity dependence of the if -band PL-relation. 
Feast (2004b) discusses the zeropoint of the Mira if -band PL- 
relation, and adopting the slope observed in the LMC (—3.47) 
derives a zeropoint of 1.00 ± 0.08, averaging over indepen- 
dently derived ZPs from trigonometric parallaxes, OH VLBI 
expansion parallaxes and Galactic Globular Clusters. Adopting 
a slope of —3.47 and refitting the ZP of the Bulge sample, the 
bias corrected value becomes 15.73 ± 0.03, and without metal- 
licity correction (consistent with the assumption above about 
the metallicities in Bulge, LMC, SMC) leads to a distance to 
the GC of 8.8 ± 0.4 kpc. This independent distance estimate is 
in between the values derived using no or a strong metallicity 
dependent zero point. 
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3 that is, M K = M K (rei) + /3\og(Z/Z re {), where M K (ref) is 
the known magnitude is a galaxy with metallicity Z Ta {, and Mk the 
magnitude it would have in a galaxy with metallicity Z 
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Appendix A: The light curve analysis model 

Some small changes w.r.t. the implementation of the lightcurve analy- 
sise model in G04 are described. 

The only change in the first part of the code is the level at which 
a period is accepted as significant. This level was set at significance 
= 1.0 x 10" 16 , compared to 5.5 x 10~ n in G04. This was possible 
as-contrary to G04-only objects with large amplitudes were searched 
for. The resulting increase in the number of spurious periods was then 
caught in the process of visual inspection. 

As in G04 a list of known objects (both known non-LPVs, and 
known AGB giants and LPVs) was compiled to ease automatic asso- 
ciation. The list comprises: 

(1) 14833 IRAS sources within 10 deg. radius of the centre of the 
49 OGLE fields at RA = 268.87, Dec = -31. 03 4 , 

(2) 51141 ISOGAL sources from Omont et al. (2003, hereafter 
OGA03; those within the extreme values of the OGLE field boundaries 
261.612 <= RA <= 276.133, -40.726 <= Dec <= -21.328), 

(3) 268 pulsating variable stars, 1650 Eclipsing Binaries, and 943 
Miscellaneous variable stars from OGLE-I (Udalski et al. 1994, 1995a, 
1995b, 1996, 1997), 

(4) 332 objects from Alard et al. (2001, hereafter ABC01) who 
correlated ISOGAL sources within the NGC 6522 and Sgr I Baade 
windows with the MACHO database, 

(5) 2353 objects from Ojha et al. (2003, hereafter OOS03) who 
studied sources in 9 ISOGAL fields, 

(6) 174 M-giants later than spectral type M0 in the NGC 6522 
Baade's window from Glass & Schultheis (2002, hereafter GS02), 

(7) 421 objects from Glass et al. (2001, and erratum), who mon- 
itored in K over four years an 24x24 arcmin area near the Galactic 
Centre, 

(8) 122 objects from Alard et al. (1996, identified as Ter- 
[number]) who identified LPV using red photographic plates, 

(9) 494 objects from Lloyd-Evans (1976, hereafter Lloyd-Evans 
(1976) who identified Mira variables in the three Baade windows, 
identified as TLE- [field] -[number] in Table |2j, Blanco et al. (1984; 
herafter BMB, who identified M giants in Baade's window), and 
Blanco (1984; herafter B84, who identified RR Lyrae variables) with 
coordinates listed in the SIMBAD database to 1" or better, 

(10) 33 Nova related objects from Cieslinski (2003). The total 
number of sources used in the automatic correlation is 72764. 

4 Retrieved from the infrared science archive at IPAC, 
http://irsa.ipac.caltech.edu/ 
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Fig. B.l. For angles of (j> of 43 and 79 degrees (left and right 
panel), the projection (drawn to scale) on the Galactic Plane of 
1300 randomly drawn stars. Disk stars are represented by filled 
circles, Bulge stars by open circles. The Sun is at X= 0, Y= 
0. The Galactic Centre is at X= 0, Y= 8.5. The lines illustrate 
Galactic longitudes of ±10 degrees. 

Appendix B: Simulations of the Galactic Bulge 
and foreground disk stars 

In this Appendix the calculations are described to model a population 
in the direction of the Galactic Bulge. 

The basic model is essentially the one proposed by Binney et al. 
(1997) to model the dust-corrected near-infrared COBE/DIRBE sur- 
face brightness map of the inner galaxy. The number density of bulge 
stars is assumed to be: 



fh = fo exp(-a /a m ) / (1 + a/a f 



(B.l) 



with f = 624, a m = 1.9 kpc, a = 0.10 kpc, f3 = 1.8. 
Binney et al. assumed a tri-axial bulge with axial ratios 1 : 0.6 : 
0.4. For numerical convenience a prolate ellipsoid is assumed here: 
a = \J x 2 + (y/ri) 2 + (z/-q) 2 with the value of 77 = 0.5 taken from 
Binney et al. 

The number of Bulge objects up to a radius r from the centre, that 
defines the probability density function in the simulation, is approxi- 
mated as: 



N b (r)= / 4na 2 f h (a)da 



(B.2) 

up to a maximum radius that is taken to be the co-rotation radius, with 
a default value of R CI = 2.4 kpc, following Dwek et al. (1995). 
The number density of disk stars is assumed to be: 

/ d = (exp(- I z I /zo) + aexp(- I z I /zt)) x 

R d (exp(-r/i? d ) - f h exp(-r/i4)) (B.3) 



with z = 0.210 kpc, z\ = 0.042 kpc, a = 0.27, R d = 2.5 kpc 
(Binney et al.) and = 1.3 kpc (Picaud & Robin 2004). This func- 
tional form follows Binney et al., but also allows for a "hole" in inner 
disk (the scaling parameter is < fh < 1, and identical to zero in 
Binney et al.). The total number of disk stars, and the probability den- 
sity functions, are defined as: 

iV d EE N d ,z(z) X iVd.rW 

given by, 



N d = 



2 / exp(— 2/20) + a exp(— 2/21) dz 
Jo 



27rr7? d (exp(-r/i?d) - fh exp(-r/R h )) dr 



(B.4) 



up to maximum values 2 max = (R CI rf), and i? max = 8.0 kpc, re- 
spectively. 

A disk or bulge star is generated according to the ratio iVb / (Nh + 
Nd,2, x Nd.r)- In case of a disk star, its height above the plane, 2, 
distance to the GC, r, and a random angle between and 2n in 
the Galactic plane, are drawn according to the probability functions 
Nd,z(z) and Nd,r{r). Its coordinates x, y, 2 are then known. 

In case of a bulge star, the distance, a, to the GC is drawn accord- 
ing to the probability function N\>, and then a star is randomly placed 
on the surface of the appropriate ellipsoid, to find x,y,z. These values 
are then rotated by an angle <j> in the Galactic plane (see Fig.|8j. 

The Galactic coordinates are then derived assuming a distance 
from the Sun to the GC of Ro = 8.5 kpc, and height above the plane 
of 20 = +24 pc (Maiz-Apellaniz 2001) 

In a second step, for every star, the known distance to the Sun is 
used to calculate its appararent magnitude, assuming an arbitrary Mk 
of —7.5 mag with a Gaussian dispersion of ok = 0.15 mag. This is 
about the dispersion observed in the PL-relation in LMC Miras when 
multi-epoch photometry is available to accurately determine mean- 
light magnitudes (Feast et al. 1989). 

In a third step, for every simulated star it is verified if it is located 
within one of the 40 lines-of-sight considered, listed in Table [6] The 
field sizes of 14.2' x 57' are approximated by a circle of radius 0.27 de- 
grees. If so, it is assumed the star would have been "detected" (Given 
the relative brightness of the LPVs it is assumed that completeness is 
not an issue). 

The number of stars drawn is such that a total of about 1200 ob- 
jects is "detected", similar to the actual number. At the end of the 
simulation, the average magnitude and dispersion per line-of-sight is 
determined, and a weighted least-square fit is made of the mean mag- 
nitude versus longitude, for all fields, and for those with | / |< 5°, as 
for the observations. In addition, the mean magnitude and dispersion 
for all "detected" stars is determined. 

Such a simulation is repeated 1000 times. Then, distribution func- 
tions and from that median and 1-sigma values of the following pa- 
rameters are determined: the number of stars (total, disk, bulge), mean 
magnitude and sigma (for every line-of-sight), mean magnitude and 
sigma for all stars, and slope and error in the slope both when fitted 
over all longitudes, and for those fields with | / |< 5°. 

For the standard model of Binney et al. described above (i.e. 
fh = 0) it turns out that for two values of a slope (because of 
the contamination by disk stars in the outer fields, the slope fitted over 
\ I I < 5° is considered for now on) in agreement with observations is 
found: cj> = 43 and 79 degrees (with values between 25 and 85 degrees 
resulting in predicted slopes within l-cr of the observed one.) 

Fieures IB. II IB. 21 and IB. 31 show the results for these two cases. 
Fi sure |B~T1 show s the distribution on the Galactic Plane for a random 
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sub-sample of all stars simulated, and illustrates a fundamental differ- 
ence between the two cases. For large viewing angles the outer fields 
/ I £10° will be dominated by disk stars. Figure lB~2l shows for the 
same random sub-sample the observed magnitude as a function of I. 
In Figure lB"31 the simulated mean magnitude and error for each field 
are compared to the observations in the top panel, while in the bottom 
panel the observed and predicted number of objects are compared. It 
is from this plot that one may conclude that the model with <j> = 43° 
is to be preferred over the one with <j> — 79° as the latter model pre- 
dicts too few stars, especially in the outer fields. Comparing only the 
observed and predicted number of stars (in a \ 2 sense) a best fit is 
found for <j> — 35° (with values between and 60 degrees resulting 
in a reduced x 2 within 1 unit of the minimum y 2 ). Combining the 
constraints from the slope and the number of stars a viewing angle of 
4> = 43 ± 17 degrees is the preferred value. 

One may consider the ratio of Bulge to Disk stars as uncertain, and 
therefore, a model was considered with <j> = 79° and fo — 350. The 
latter value was set so that the model predicted the observed number 
of stars in the I — —11° field. Such a model would still underestimate 
the number of stars in the outer fields at positive I, and would give a 
slope no longer in agreement with observations. 

Finally, a model including a hole in the inner disk is considered 
(i.e. /h = 1). To have the same ratio of bulge to disk stars, fo was set 
to 425. The results are very similar and the best fitting angle is now 40 
degrees. 

For reference, the predicted number of disk and bulge stars for the 
viewing angles of 43 and 79 degrees, and for the model with cj> = 40° 
and the central hole in the disk, are listed in Tab. IB. II As they are 
quite different, these predictions may be usefull when additional data 
(proper motion 5 , radial velocities) become available to constrain the 
ratio of disk to bulge objects as a function of galactic coordinates. 

Appendix C: Comparing stellar evolution codes 

From P. Wood's webpage l http://www.mso.anu.edu.au/~wood/ 1 the 
models described in Vassiliadis & Wood (1993, VW) were down- 
loaded. These files list for the individual calculated models on the 
AGB the relevant stellar quantities (remaining stellar mass, luminosity 
and effective temperature amongst other quantities) and the evolution- 
ary time. They are available for Z = 0.016 (1.0, 1.5, 2.0, 2.5, 3.5, 5.0 
Mq), Z = 0.008 (0.945, 1.0, 1.5, 2.0, 2.5, 3.5, 5.0 M ), Z = 0.004 
(0.89, 1.0, 1.5, 2.0, 2.5, 3.5, 5.0 M ) and Z = 0.001 (1.0, 1.5 M ). 

For our comparison with simulation we used the solar metallic- 
ity models, which are expected to be the most representative for our 
Bulge sample. However, different studies indicate that the Bulge may 
have quite a broad metallicity distribution, peaking at about —0.25 dex 
with dispersion of 0.3 dex (see e.g. McWilliam & Rich 1994, Ramirez 
et al. 2000). The AGB lifetime, LPV lifetime and LPV period distri- 
bution was determined. The AGB lifetime is defined as the time be- 
tween the first model in the file (the start of the AGB) upto the point 
where the remaining envelope mass becomes less than 0.04 Mq, or 
left > 4500A", that is taken as the start of the post- AGB evolution. 

For each timestep the fundamental period is calculated following 
VW. The star is assumed to be in the Mira instability strip when the 
bolometric magnitude is within 0.20 magnitude (the assumed width of 
the instability strip at a given period) of the PL-relation (Feast et al. 
1989): 

M bol = -3.00 log P + 2.85 (C.l) 

5 Sumi et al. (2004) present proper motions for 5 million objects in 
the OGLE fields, centered on the expected Jo magnitude and (V — I)o 
colour of red clump giants but also including some red giants. Of the 
2691 LPVs in the present study, 1612 are listed in Sumi et al. 



assuming a LMC distance of 18.50, and when the mass loss rate is be- 
low a critical value, as this PL-relation is derived for optically visible 
objects and the LPV samples studied in G04 and in the present paper 
have been culled by only considering objects with (J — K)o < 2.0. In 
such a way the lifetime and the period distribution of optically visible 
LPVs can be determined. 

The critical mass loss rate is determined by taking for each mass 
in the grid of solar metallicities typical values of luminosity and effec- 
tive temperature inside the instability strip and then using a radiative 
transfer program (Groenewegen 1993) with the appropriate model at- 
mosphere for M-stars (Fluks et al. 1994), and typical dust properties 
(silicate dust, condensation temperature of 1500 K, dust-to-gas ratio 
0.005, expansion velocity 10 km s _1 ) to determine the critical mass 
loss rate at which the star would become redder in (J — K) than 
2.0. The critical mass loss rates found are between 4 and 20 x 10 -6 
M yr _1 depending on the initial mass of the model. In fact, the criti- 
cal mass loss rate is observed to scale with \/~L, as expected as the dust 
optical depth predominantly determines the infrared colours, hence the 
mass loss rate is proportional to the stellar radius, all other things be- 
ing equal. For the 1.0 and 1.5M initial mass models the mass loss 
rate inside the instability strip always remains below the critical mass 
loss rate. The adopted critical mass loss rate is 1.0 xlO" 5 ^L/lWm 
M yr _1 . The results for the solar metallicty models are summarised 
in TablelcTl 

We will now consider the synthetic AGB models of Wagenhuber 
& Groenewegen (1998; WG). The reason being that the VW models 
exist only for a limited number of initial masses and secondly be- 
cause mass loss on the RGB was only included for initial masses be- 
low 1M , while it is well known that the effect of mass loss on the 
RGB is substantial also at and above 1 M . The VW mass loss rate 
recipe was implemented, and to mimic the VW tracks as closely as 
possible the mixing length parameter a (basically setting the effective 
temperature scale) in the WG models was tuned to give similar AGB 
and LPV lifetimes. The results are summarised in Table lC. fl It turns 
out that with a — 2.0 the lifetimes match very well especially at low 
initial mass. 

Mass loss on the RGB is described by a Reimers law with a scal- 
ing factor i)rgb = 0.35. This gives the required mass loss (0. 13, 0. 16, 
0.17 M for a star of l.OM initial mass at Z = 0.004, 0.008, 0.019, 
respectively; M. Salaris, private communication) to give the observed 
mean colour of Horizontal Branch stars in Galactic Globular Clusters. 
Table ?? summarises for a set of WG models with ?)rgb = 0.35 and 
a — 2.0 the AGB and LPV lifetime, and Figure IcTl displays the pe- 
riod distribution of optically visible stars inside the instability strip 
(normalised to one each time) for a few initial masses. 

Table ?? also includes for a few initial masses the results if 
slightly different metallicities of Z= 0.01 and 0.02 are adopted, and 
Figure lG2l shows the corresponding Period distribution. These results 
indicate that the effect of metallicity on the pulsation properties for 
the typical metallicities in the Bulge is small. 
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Fig. B.2. For angles of of 43 and 79 degrees (left and right panel), the distribution of A'-magnitudes as a function of Galactic 
longitude for the stars shown in Fig. IB. II Disk stars are represented by filled circles, Bulge stars by open circles. 
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Fig. B.3. For angles of <j> of 43 and 79 degrees (left and right panel), the top panel compares the observed (open squares) and 
modelled (filled circles) mean magnitudes and errors for the 40 fields listed Table[6] with the line indicating the best fit over the 
range | I \ < 5°. For clarity the observed and modelled points are slightly off-setted in longitude. The bottom panel compares the 
observed (open squares) and modelled (filled circles) number of stars per field. For both angles the observed slope over the range 
| / |< 5° is reproduced, but the distribution of stars is much better fitted for cf> = 43°. 
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Table 4. Properties of the OGLE-fields 



BUL_SC 


I 


b 


Total (o) 


Variable (b) 




LPV (d) 


1 


1.08 


-3.62 


730 


4597 


1.68 / 1.49 


42 


2 


2.23 


-3.46 


803 


5279 


1.55 / 1.65 


48 


3 


0.11 


-1.93 


806 


8393 


2.89 


115 


4 


0.43 


-2.01 


774 


9096 


2.59 / 2.94 


86 


5 


-0.23 


-1.33 


434 


7257 


5.73 /-/4.13 


118 


6 


-0.25 


-5.70 


514 


3211 


1.37 


47 


7 


-0.14 


-5.91 


463 


1618 


1.33 / 1.28 


21 


8 


10.48 


-3.78 


402 


2331 


2.14 


8 


9 


10.59 


-3.98 


330 


1847 


2.08 


21 


10 


9.64 


-3.44 


458 


2499 


2.23 


16 


1 1 


9.74 


-3.64 


426 


2256 


2.27 


18 


12 


7.80 


-3.37 


535 


3476 


2.29 / 2.20 


33 


13 


7.91 


-3.58 


570 


3084 


2.06/1.82 


21 


14 


5.23 


2.81 


619 


4051 


2.49 


51 


15 


5.38 


2.63 


601 


3853 


2.77 


71 


16 


5.10 


-3.29 


700 


4802 


2.15 /2.23 


45 


17 


5.28 


-3.45 


687 


4690 


1.94/2.29 


167 


18 


3.97 


-3.14 


749 


5805 


1.83 


55 


19 


4.08 


-3.35 


732 


5255 


2.74 


51 


20 


1.68 


-2.47 


785 


5910 


1.94/2.02 


64 


21 


1.80 


-2.66 


883 


7449 


1.83 / 1.78 


60 


22 


-0.26 


-2.95 


715 


5589 


2.74 


70 


23 


-0.50 


-3.36 


723 


4815 


2.70 


60 


24 


-2.44 


-3.36 


612 


4304 


2.52 


56 


25 


-2.32 


-3.56 


622 


3046 


2.34 


61 


26 


-4.90 


-3.37 


728 


4713 


1.86 


39 


27 


-4.92 


-3.65 


691 


3691 


1.69 


46 


28 


-6.76 


-4.43 


406 


1472 


1.64 


16 


29 


-6.64 


-4.62 


492 


2398 


1.53 


36 


30 


1.94 


-2.84 


762 


6893 


1.91 / 1.78 


50 


31 


2.23 


-2.94 


790 


4789 


1.81 / 1.74 


59 


32 


2.34 


-3.14 


797 


5007 


1.61 / 1.82 


54 


33 


2.35 


-3.66 


739 


4590 


1.70 / 1.82 


105 


34 


1.35 


-2.40 


961 


7953 


2.52/2.32 


60 


35 


3.05 


-3.00 


771 


5169 


1.84/2.20 


39 


36 


3.16 


-3.20 


873 


8805 


1.62/1.52 


38 


37 


0.00 


-1.74 


664 


8367 


3.77 


108 


38 


0.97 


-3.42 


710 


5072 


1.83 / 1.94 


57 


39 


0.53 


-2.21 


784 


7338 


2.63 / 2.70 


99 


40 


-2.99 


-3.14 


631 


4079 


2.94 


55 


41 


-2.78 


-3.27 


603 


4035 


2.65 


49 


42 


4.48 


-3.38 


601 


4360 


2.29 


32 


43 


0.37 


2.95 


474 


3351 


3.67 


112 


44 


-0.43 


-1.19 


319 


7836 


6.0/-/ 6.00 


132 


45 


0.98 


-3.94 


627 


2262 


1.64/1.53 


32 


46 


1.09 


-4.14 


552 


2057 


1.71/1.65 


26 


47 


-11.19 


-2.60 


301 


1152 


2.60 


12 


48 


-11.07 


-2.78 


287 


973 


2.35 


12 


49 


-11.36 


-3.25 


251 


826 


2.09 


18 


total 






30490 


221701 




2691 



(a) Total number of objects detected in the field. From Udalski et al. (2002), in units of 10 3 objects 

(b) Total number of candidate variable stars. From Wozniak et al. (2002). 

(c) Visual extinction. From Sumi (2004), except for SC44, where Ay = 6.0 has been adopted based on the proximity to SC5. The second 
value-when listed-comes from Popowski et al. (2003). The third value-when listed-comes from Schultheis et al. (1999b). 

(d) Total number of LPVs. 
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Table 5. Current Period compared to values listed in the literature 



OGLE-name 



Periods 



Remark 



harmonic ? 
harmonic ? 



harmonic ? 



bul_sc01_0558 232.2 (OGLE-II), 217.0, Mira (OGLE-I) 

bul.sc01_1467 519.3 (OGLE-II), 265.0, Mira (OGLE-I) 

bul.sc01.1738 458.0 (OGLE-II), 227.5, Mira (OGLE-I) 

bul.sc01_2079 524.7 (OGLE-II), 520.0, Mira (OGLE-I) 

bul.sc20_0832 245.9 (OGLE-II), 254.1: (ABC01), 235 (LE76), 336 (GWCF) 

bul.sc20.0975 167.7 (OGLE-II), 145.8 (ABC01) 

bul.sc20_l 189 220.0 (OGLE-II), 121.3 (ABC01) 

bul.sc20_1292 87.7 (OGLE-II), 88.1 (ABC01) 

bul.sc20_1761 231.6 (OGLE-II), 107.9 (ABC01), 240 (LE76), 235 (GWCF) 

bul.sc20_l 826 33 1 .7 (OGLE-II), 330: (LE76) 

bul.sc20_1928 297.6 (OGLE-II), 300 (LE76), 293 (GWCF) 

bul.sc20_2013 453.3 (OGLE-II), 474.2: (ABC01), 500:: (LE76), 480 (GWCF) 

bul.sc20.2269 409.3 (OGLE-II), 400: (LE76), 383 (GWCF) 

bul.sc20.2291 311.2 (OGLE-II), 306.9: (ABC01), 315 (LE76), 308 (GWCF) 

bul.sc34.3759 289.8 (OGLE-II), 293.8 (ABC01), 265 (LE76), 237 (GWCF) 

bul.sc45.0704 405.3 (OGLE-II), 430 (LE76), 467, Mira (OGLE-I), 413 (GS03) 

bul.sc45.1068 115.2 (OGLE-II), 115 (LE76), 116 (GCVS), 116 (GS03) 

bul.sc45.1586 193.3 (OGLE-II), 190 (LE76), 195 (GS03) 

bul.sc46.0866 322. 1 (OGLE-II), 305 (LE76), 320 (GS03) 

bul.sc46_l 163 343.7 (OGLE-II), 329.0, Mira (OGLE-I) 
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Table 6. Zero point of the LT-band PL-relation. 



Fields 


jl<0 


6 l«) 




N (c) 




(deg) 


(deg) 


(mag) 




8,9,10,11 


9.98 


-3.75 


15.58 ± 0.59 


23 


12,13 


7.74 


-3.50 


15.43 ± 0.44 


27 


15 


5.34 


+2.59 


15.38 ± 0.48 


22 


17 


5.21 


-3.50 


15.36 ± 0.47 


24 


14 


5.14 


+2.73 


15.36 ± 0.38 


29 


16 


5.04 


-3.37 


15.41 ± 0.62 


17 


42 


4.29 


-3.52 


15.46 ± 0.45 


13 


19 


3.98 


-3.41 


15.35 ± 0.28 


20 


18 


3.90 


-3.19 


15.37 ± 0.45 


20 


35,36 


3.04 


-3.18 


15.42 ± 0.39 


29 


33 


2.24 


-3.72 


15.55 ± 0.37 


22 


32 


2.26 


-3.22 


15.57 ± 0.37 


24 


31 


2.15 


-3.01 


15.52 ± 0.41 


32 


2 


2.08 


-3.59 


15.49 ± 0.39 


19 


30 


1.90 


-2.90 


15.56 ± 0.41 


25 


21 


1.77 


-2.70 


15.50 ± 0.44 


33 


20 


1.64 


-2.52 


15.42 ± 0.43 


41 


34 


1.33 


-2.43 


15.29 ± 0.33 


35 


46 


1.08 


-4.14 


15.46 ± 0.32 


18 


1 


0.94 


-3.69 


15.51 ± 0.34 


22 


45 


1.00 


-3.96 


15.44 ± 0.51 


15 


38 


0.88 


-3.49 


15.44 ± 0.36 


25 


39 


0.48 


-2.26 


15.30 ± 0.39 


69 


4 


0.30 


-2.10 


15.48 ± 0.46 


54 


43 


0.21 


+2.84 


15.42 ± 0.35 


47 


3 


-0.01 


-2.02 


15.43 ± 0.40 


70 


37 


-0.09 


-1.81 


15.34 ± 0.38 


68 


7 


-0.09 


-5.86 


15.67 ± 0.36 


6 


5 


-0.31 


-1.39 


15.38 ± 0.42 


90 


6 


-0.28 


-5.70 


15.80 ± 0.46 


8 


22 


-0.34 


-3.01 


15.50 ± 0.42 


46 


44 


-0.48 


-1.26 


15.36 ± 0.36 


79 


23 


-0.63 


-3.44 


15.53 ± 0.40 


30 


25 


-2.41 


-3.64 


15.44 + 0.31 


33 


24 


-2.56 


-3.46 


15.62 ± 0.34 


29 


41 


-2.86 


-3.35 


15.63 ± 0.38 


24 


40 


-3.16 


-3.27 


15.50 ± 0.37 


34 


26,27 


-4.99 


-3.58 


15.63 ± 0.42 


39 


28,29 


-6.75 


-4.55 


15.55 ± 0.42 


13 


47,48,49 


-11.19 


-2.78 


15.95 ± 0.35 


18 



(a) Mean values of individual objects 

(b) Zero point of the P-band PL-relation for a slope of —3.37. 

(c) Number of objects in Box "C" used to calculate the PL-relation. 
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Table B.l. Observed total, and predicted number of disk and bulge stars per field. 









d> = 


43° 


d> = 


79° 


= 4O°,/ h = 1 


1 


b 


Nobs 


-/Vdisk 


iVbulgc 


-/Vdisk 


TVbulgc 


-/Vdisk 


TVbulgc 


9.98 


-3.75 


5.75 


2.1 ± 1.5 


3.1 ± 2.0 


2.8 ± 0.9 


0.0 ± 0.0 


2.0 ± 1.5 


3.9 ± 1.9 


7.74 


-3.50 


13.5 


3.1 ± 1.6 


6.8 ± 2.6 


2.8 ± 1.6 


1.8 ± 1.2 


3.0 ± 1.7 


7.5 ± 2.8 


5.34 


+2.59 


22 


5.6 ± 2.4 


14.0 ± 3.5 


4.0 ± 1.9 


7.4 ± 2.6 


4.2 ± 2.2 


15.4 ±3.9 


5.21 


-3.50 


24 


3.8 ± 1.7 


12.0 ± 3.6 


2.9 ± 1.7 


6.6 ± 2.5 


2.9 ± 1.7 


12.9 ± 3.6 


5.14 


+2.73 


29 


5.0 ± 2.3 


13.9 ± 3.4 


3.8 ± 1.8 


7.9 ± 2.8 


3.9 ± 2.0 


15.3 ±4.0 


5.04 


-3.37 


17 


3.9 ± 2.0 


13.0 ± 3.5 


2.9 ± 1.7 


7.4 ± 2.7 


3.0 ± 2.0 


13.9 ± 3.5 


4.29 


-3.52 


13 


3.8 ± 1.8 


14.0 ± 3.6 


2.8 ± 1.7 


9.5 ± 3.0 


2.9 ± 1.5 


15.2 ± 3.7 


3.98 


-3.41 


20 


3.9 ± 2.0 


15.5 ± 4.1 


2.9 ± 1.8 


11.2 ± 3.4 


3.0 ± 1.8 


16.6 ±4.0 


3.90 


-3.19 


20 


4.8 ± 2.2 


17.2 ± 4.0 


3.7 ± 1.7 


12.7 ± 3.4 


3.7 ± 1.8 


18.7 ±4.1 


3.04 


-3.18 


14.5 


4.9 ± 2.2 


19.9 ± 4.4 


3.8 ± 1.7 


17.1 ± 4.1 


3.8 ± 1.7 


21.1 ±4.7 


2.24 


-3.72 


22 


3.0 ± 1.9 


16.9 ± 4.1 


2.1 ± 1.5 


17.4 ± 4.1 


2.0 ± 1.5 


18.3 ± 3.9 


2.26 


-3.22 


24 


4.9 ± 2.3 


21.7 ± 4.5 


3.8 ± 2.0 


21.8 ± 4.4 


3.1 ± 1.9 


23.5 ±4.7 


2.15 


-3.01 


32 


5.7 ± 2.2 


24.9 ± 5.0 


4.0 ± 1.9 


24.5 ± 4.8 


3.9 ± 1.8 


26.5 ± 4.8 


2.08 


-3.59 


19 


3.9 ± 1.9 


18.5 ± 4.3 


2.9 ± 1.7 


19.4 ± 4.3 


2.8 ± 1.5 


19.9 ±4.5 


1.90 


-2.90 


25 


6.1 ± 2.5 


26.7 ± 5.0 


4.7 ± 2.1 


27.9 ± 5.1 


4.1 ± 2.0 


28.7 ± 5.0 


1.77 


-2.70 


33 


7.1 ± 2.6 


30.4 ± 5.4 


5.2 ± 2.1 


31.6 ± 5.4 


5.0 ± 2.3 


32.2 ± 5.4 


1.64 


-2.52 


41 


8.3 ± 2.9 


33.8 ± 5.6 


6.0 ± 2.4 


36.1 ± 6.0 


5.8 ± 2.2 


35.9 ± 6.0 


1.33 


-2.43 


35 


8.8 ± 2.8 


36.6 ± 6.1 


6.6 ± 2.5 


41.0 ± 6.3 


5.9 ±2.5 


39.1 ±5.9 


1.08 


-4.14 


18 


2.2 ± 1.5 


15.3 ± 3.8 


1.9 ± 1.3 


17.4 ± 4.0 


1.9 ± 1.3 


16.5 ± 3.9 


0.94 


-3.69 


22 


3.7 ± 1.9 


19.6 ± 4.4 


2.8 ± 1.6 


22.9 ± 4.6 


2.0 ± 1.5 


21.0 ±4.2 


1.00 


-3.96 


15 


2.9 ± 1.7 


17.0 ± 4.0 


2.0 ± 1.4 


19.6 ± 4.3 


2.0 ± 1.3 


18.1 ±4.1 


0.88 


-3.49 


25 


4.0 ± 1.9 


21.7 ± 4.6 


3.0 ± 1.7 


25.4 ± 4.8 


2.9 ± 1.7 


23.4 ± 4.5 


0.48 


-2.26 


69 


10.1 ± 3.0 


43.3 ± 6.4 


7.5 ± 2.6 


53.2 ± 7.2 


6.9 ± 2.6 


46.0 ± 6.7 


0.30 


-2.10 


54 


11.6 ± 3.1 


48.1 ± 6.3 


8.6 ± 2.8 


59.6 ± 7.2 


7.9 ± 2.7 


51.0 ±6.6 


0.21 


+2.84 


47 


5.0 ± 2.4 


26.9 ± 5.0 


3.9 ± 2.1 


31.8 ± 5.2 


3.7 ± 1.7 


28.2 ± 5.1 


-0.01 


-2.02 


70 


12.0 ± 3.5 


50.6 ± 6.4 


9.0 ± 2.9 


63.6 ± 7.0 


8.5 ± 2.7 


53.6 ±6.8 


-0.09 


-1.81 


68 


14.2 ± 3.6 


57.8 ± 7.4 


10.7 ± 3.0 


73.3 ± 8.1 


10.1 ± 3.1 


61.6 ±7.5 


-0.09 


-5.86 


6 


1.0 ± 0.8 


6.0 ± 2.4 


0.9 ± 0.7 


7.0 ± 2.5 


0.0 ± 0.5 


6.3 ± 2.6 


-0.31 


-1.39 


90 


20.7 ± 4.4 


77.4 ± 8.5 


15.2 ± 3.8 


98.4 ± 9.7 


14.8 ± 3.8 


81.9 ± 8.6 


-0.28 


-5.70 


8 


1.0 ± 0.8 


6.9 ± 2.6 


0.9 ± 0.7 


7.7 ± 2.8 


0.0 ± 0.7 


7.1 ± 2.6 


-0.34 


-3.01 


46 


5.8 ± 2.2 


29.1 ± 5.4 


4.2 ± 2.0 


35.0 ± 5.5 


3.9 ± 1.8 


30.4 ± 5.6 


-0.48 


-1.26 


79 


23.3 ± 4.7 


83.5 ± 8.8 


17.4 ± 4.1 


105.5 ± 9.9 


16.7 ± 4.0 


88.6 ± 9.0 


-0.63 


-3.44 


30 


4.0 ± 1.9 


22.9 ± 4.5 


3.0 ± 1.8 


27.1 ± 5.0 


2.9 ± 1.6 


24.4 ± 4.9 


-2.41 


-3.64 


33 


3.1 ± 1.7 


18.3 ±4.2 


2.1 ± 1.5 


17.5 ±3.9 


2.1 ± 1.5 


19.8 ±4.4 


-2.56 


-3.46 


29 


3.9 + 2.1 


19.9 ±4.1 


2.9 ± 1.7 


17.9 ±4.1 


2.9 ± 1.7 


21.2 ±4.6 


-2.86 


-3.35 


24 


4.1 ±2.1 


20.1 ±4.2 


3.0 ± 1.9 


17.1 ±4.2 


3.0 ± 1.8 


21.2 ±4.5 


-3.16 


-3.27 


34 


4.3 ± 2.0 


19.8 ±4.3 


3.1 ± 1.7 


16.1 ±4.1 


3.0 ± 1.8 


21.4 ±4.5 


-4.99 


-3.58 


19.5 


3.1 ± 1.7 


12.2 ± 3.5 


2.8 ± 1.6 


6.1 ± 2.6 


2.9 ± 1.6 


13.5 ±3.5 


-6.75 


-4.55 


6.5 


1.9 ± 1.2 


4.9 ± 2.2 


2.0 ± 1.0 


0.9 ± 0.7 


1.0 ± 1.0 


5.2 ±2.1 


-11.19 


-2.78 


6.0 


4.9 ± 2.0 


0.0 ± 0.0 


3.9 ± 1.5 


0.0 ± 0.0 


5.0 ± 2.0 


0.0 ± 0.5 
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Table C.l. AGB lifetime and optically visible LPV lifetimes for two sets of models 



Vassiliadis & Wood 



Wagenhuber & Groenewegen 



z 


Mass 


MS lifetime 


AGB lifetime 


LPV lifetime 


AGB lifetime 


LPV lifetime 


remark 




(M©) 


(10 9 years) 


(10 3 years) 


(10 3 years) 


(10 3 years) 


(10 3 years) 




0.016 


1.0 


11.25 


595 


101 


487 


49 


JlRPB = 0. O-— 1.9 












560 


93 


r? R pi R — ry— ? 












595 


129 


?]R,GB — 0, Oi— 2.1 


0.016 


1.5 


2.742 


929 


272 


873 


303 


71t) r< r — f) Cv— 1 Q 












942 


284 


77 RPR = 0, a— 2.0 












1019 


282 


riDpn — Q CV— ? 1 


0.016 


2.0 


1.236 


1168 


153 
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Fig. C.l. Theoretical period distribution of optically visible stars inside the observed instability strip for masses 1.1, 1.2, 1. 
2.0, 2.5, 3.0M o (left to right, top to bottom). 
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Fig. C.2. As Figure ICTl illustrating the influence of metallicity on the theoretical distribution. In the left hand panels for Z 
0.010 and on the right for Z = 0.020, for stars with inital masses 1.2 (top), 2.0 (middle) and 3.0 (bottom) M . 



